namespace hms
{

/** @brief Returns a conditional type matrix which evaluates to true for all
 * cells whose water depth is below SWESolver::minWaterDepth. */
template<typename MeshType>
auto SWESolver<MeshType>::isDry() const -> decltype(auto) {
	return d().values() < minWaterDepth;
}

/** @brief Returns a conditional type matrix which evaluates to true for all
 * domain cells whose water depth is below SWESolver::minWaterDepth */
template<typename MeshType>
auto SWESolver<MeshType>::isDryDomain() const -> decltype(auto) {
	return d().domainValues() < minWaterDepth;
}

} // end namespace hms
